Theory of Nonlinear Magneto-Optics 



U. Pustogowa, T. A. Luce, W. Hiibner, and K. H. Bennemann 
Institute for Theoretical Physics, Freie Universitdt Berlin, Arnimallee I4, D - 14195 Berlin, 

Germany 

Abstract 

Nonlinear magneto-optics is a very sensitive fingerprint of the electronic, 
magnetic, and atomic structure of surfaces, interfaces, and thin ferromagnetic 
films. Analyzing theoretically the nonlinear magneto-optical Kerr effect for 
thin films of Fe(OOl) and at Fe surfaces we demonstrate exemplarily how vari- 
ous electronic material properties of ferromagnets , such as the d-band width, 
the magnetization, the substrate lattice constant, and the film-thickness de- 
pendence can be extracted from the calculated nonlinear Kerr spectra. Fur- 
thermore, we show how the substrate d electrons (Cu(OOl)) affect the nonlin- 
ear Kerr spectra even without being themselves spin-polarized and without 
film-substrate hybridization. We show that the Kerr rotation angle in second 
harmonic generation is enhanced by one to two orders of magnitude compared 
to the linear Kerr angle and how symmetry can be used to obtain the direc- 
tion of magnetization in thin films and at buried interfaces from nonlinear 
magneto-optics. 
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I. INTRODUCTION 



Nonlinear magneto-optics for metals is a new expanding field the significance of 
which results from extending electrodynamics to the nonlinear case and from applications. 
Particularly, an important task of theory was to extend the symmetry analysis to an elec- 
tronic theory since this allowed second harmonic generation (SHG) to be used as an 
important new tool for the analysis of surface [§, interface, thin film, and multilayer 
magnetism. The surface sensitivity results since breakdown of inversion symmetry is nec- 
essary for SHG. Thus, it is possible to determine with the help of SHG surface magnetic 
moments, magnetic anisotropy, domains, magnetostriction, and, in particular, properties of 
buried interfaces. 

Using an electronic theory (with a bandstructure input) we calculate (i) the material spe- 
cific nonlinear response function M, As.o., c^), depending on magnetization M, spin- 
orbit coupling Xs o. and film thickness d, (ii) the nonlinear Fresnel formula for the reflected 

(2) 

nonlinear fields £^(20;), and thus finally, (iii) the complex nonlinear Kerr rotation angle ■ 
The results, presented in the following demonstrate, in particular, what can be learned (a) 
from an analysis of the frequency dependence and (b) from the symmetry analysis and the 
Kerr rotation. Some selected general results obtained previously are shown, but also inter- 
esting new results indicating important substrate effects. Remarkably, the calculation |^ 
and recent experiments |^JT^ yield that is larger by two orders of magnitude compared 



to the linear Kerr rotation. 

II. THEORY 
A. Nonlinear susceptibility 

Using response theory the nonlinear optical polarization is given by 

P^ = x'^IEj{co)Ei{u), (1) 

(2) 

where Pi is the i-th component of the polarizability and where the susceptibility tensor xl/i 



is determined by p|,p 
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Here E'kio- are the electronic energy levels resulting from the bandstructure calculations, 
(kZa"|2|k + qy, I'a) are the transition matrix elements, which control the symmetry and thus 
the interface sensitivity. If the matrix elements are treated as constants one has to introduce 
a cut-off to guarantee the surface sensitivity. This is done by the factor C [|1^]. Typically, 
magnetic dipole effects are negligible if the electric dipole contribution does not vanish. 
Xs.o. is the key parameter for nonlinear magneto-optics. In particular, for determing the 
Kerr effect it might be useful to decompose x''^^ iii odd and even contributions x^"^^ = 
^(2)(M)+x(2)(M2). 

For the application of the general theory thin films are of particular interest. Then, the 
layer-dependent contributions to x^"^^ reveal interesting film properties and substrate effects. 
For the thin film calculations we use C = Wk+2q||,/''o-Wk+q||,ro-M/kio-, where Wa denotes the 
weight of the density of state |k/cr) in the Wigner-Seitz cell of the first monolayer. Note, 
due to the factor C SHG will always result from the surface layer or from the layer at 
the interface, however the electronic structure of this layer depends on the film thickness. 
The bandstructure E'kicr is calculated employing the full-potential linear muffin-tin orbital 
method (FP-LMTO) [jT^ for thin Fe(OOl) films within a symmetric slab calculation. 

Results for one, three, five, and seven layers of Fe are presented in Fig. H The depth 
of the first minimum in the spectra scales linearly with the magnetic moment of the top 
layer while the energy of this minimum reflects the c?-band width. As can be seen from 
Fig. 1^, the nonlinear spectra reveal the structural changes induced by different substrates or 
temperatures. In Fig. ^ we compare the ab initio nonlinear Kerr spectrum for the Fe(OOl) 
monolayer with tight-binding calculations. From the results of Fig. ^ it becomes obvious 
that it is necessary to use ab initio calculations for the films if only bulk tight-binding 
parameters are available. However, tight-binding calculations may be adequate if one uses 
as input parameters those determined from ab initio calculations for thin films. 
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We expect generally also interesting results for x^^^ at heterogeneous interfaces like 
Fe/Cu(001). For example, it is of interest to find out how the substrate interferes with 
the electronic structure of the thin film. For this reason we performed calculations of the 
Fe/Cu bilayer. Results are shown in Fig. The additional peak structure reveals the elec- 
tronic influence even of the nonmagnetic Cu on the nonlinear magneto-optic response of thin 
films. This might be of general significance also with respect to the analysis of quantum 
well states (tuning of interface properties, Xs,o. in substrate, M in films). 

From the results of Figs. 1 through 4 we learn that the nonlinear Kerr effect is a sensitive 
probe of the magnetic and electronic structure at interfaces. 



B. Symmetry and Kerr rotation 



That SHG is very symmetry sensitive is furthermore particular demonstrated by its 
polarization dependence and by the remarkable enhancement of the Kerr rotation The 
latter important result was first theoretically derived by Pustogowa et al. and later 



observed in experiments by Bohmer et al. and Koopmans et al. ||T0|. 

First, we note that due to the symmetry sensitivity already SHG exhibits a correspond- 
ing polarization dependence, which interestingly is very different for transition and noble 
metals . The latter results, since s and d electrons feel the breakdown of inversion symme- 
try at the interface differently. This symmetry dependence mainly results from the matrix 
elements in x^'^\ see Eq. (0). The different possible optical excitations contributing to x*'^'* 
for noble and transition metals are illustrated in Fig. ^. Obviously, SHG yields interesting 
information on the electronic structure even for M=0 and when no Kerr effect is present. 
This area needs further analysis and might offer interesting new results Q. 

Extending the symmetry analysis to the magnetic case we determine the Kerr rotations 
for different configurations and polarizations. The general analysis can be performed by 
using for s and p polarization (see |Q) 
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where in the longitudinal configuration the susceptibility tensor is given by 
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and in the polar configuration by 
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Here, $ and ip denote the angles of polarization of the reflected frequency doubled and of 
the incident light, Ap^g are the amplitudes, fc,s and Fc^s are the Fresnel coefficients, and ts,p 
are the linear transmission coefficients for the fundamental light. Then, from 



(2) ^ 
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we determine the Kerr rotation for the polar and longitudinal configuration and for s and p 
polarized incident light. The following results are obtained: 

For the longitudinal configuration (My) and p (s) polarized incident light one gets 
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For the polar configuration {M±) and p (s) polarized incident light one gets 

= Re- 
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The coefficients ai . . . contain explicitly the whole information about the Kerr geometry 
(directions of light incidence and polarization, crystal magnetization) and about the linear 
transmission and reflection (see IQ). 

Summarizing the results, we find that generally (j)^^^^ > 0^''*^" and 0^^^, > (p^lp^ which is 



in agreement with the experiment by Koopmans et al. . Numerical results of our theory 
using as input parameters the phase ratio of the various different tensor elements from Ni and 
linear complex indices of refraction are presented in Fig. ^ These results demonstrate clearly 
that the nonlinear Kerr angle depends sensitively on the direction of the magnetization and 
the incident beam, as well as on the light polarization. Depending on the experimental 
conditions the nonlinear Kerr angle might become even as large as 90°. Moreover, also the 
ellipticity will yield correspondingly large magnetic effects. From these results we conclude 
that the enhanced nonlinear Kerr rotation allows for the determination of easy axis and 
canted-spin configurations due to large magnetic contrast of different interfaces. 

This selective discussion demonstrates the usefulness of nonlinear magneto-optics as a 
promising spectroscopy. Our results show that SHG is a new sensitive tool for determining 
magnetic interface properties such as the interplay of structure and magnetism, magnetic 
anisotropy, and magnetic reorientation. In contrast to linear optics, which probes character- 
istic film-averaged features, the nonlinear Kerr effect originates essentially from the surface 
and interfaces and thus allows the analysis of buried interfaces. 

For future studies on the nonlinear magneto-optical Kerr effect it will be interesting to 
analyze magnetic anisotropy and easy axis, magnetic structure and domains, multilayers, and 
magnetostriction effects. Nonlinear optics at interfaces is a bulk-background-free technique 
and therefore, should yield a pronounced visibility of quantum well states and enhanced 
nonlinear paramagnetic Kerr oscillations in the presence of an external field. 

An important future development of the nonlinear Kerr effect might result from femto- 
second dynamics which is able to follow the charge and spin dynamics of metallic interfaces 
in real time without involving the lattice. 
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FIGURES 

FIG. 1. Ab initio calculated nonlinear magneto-optical Kerr spectra of Fe using Eq. (2) 
for a monolayer (dashed curve - 1), and films with 3 layers (dashed-dotted curve - 2), 5 layers 
(long-dashed - 3), and 7 atomic layers (dotted curve - 4). The second harmonic response results 
from the first atomic layer. This SH response is obtained by projecting the wave functions to the 
first atomic layer yielding the factor C. 

FIG. 2. Film-lattice-constant dependence of ab initio calculated nonlinear Kerr spectra of a Fe 
monolayer. The solid curve refers to the bulk bcc Fe lattice constant a=2.76 A, the dashed curve 
to a=2.776 A (bulk Au), and the dotted curve to a=2.783 A (bulk Ag). The long-dashed curve 
refers to the experimental a=2.88 A for Fe. Choosing lattice constants a refering to Au and Ag 
should simulate substrate effects. The inset shows for an enhanced scale the effects of different 
lattice constants for the zero of Im x^^^ aX hco ^ 3 eV and indicate structural effects which should 
be observable. 

FIG. 3. Comparison of ab initio (solid line - 2) and semi-empirical calculations of the nonUnear 
magneto-optical susceptibilities for a Fe monolayer. Semi-empirical calculations are performed for 
a two-dimensional atomic configuration (dashed-dotted curve - 1) and for a three-dimensional one 
with reduced hopping parameters (dashed curve - 3) containing the reduced coordination number 
of a monolayer compared to the bulk. 

FIG. 4. Nonlinear magneto-optical Kerr-spectrum of a Fe/Cu bilayer system without (solid 
curve) and with (dashed curve) hybridization between the Fe and Cu layers. Note, the spectrum 
look rather different from the corresponding one for freestanding monolayers indicating a strong 
effect of Cu on the optical transitions. 

FIG. 5. Illustration of second harmonic generation from noble metals [(a) and (b)] and tran- 
sition metals [(c) and (d)]. For noble metals, in case (a) no d electrons can be optically excited, 
which in contrast is possible in case (b) . In case (c) for transition metals, predominantly d electrons 
contribute to the SHG yield, whereas in case (d) the excitation starts from the s band. Note, the 
cases (a) and (c) refer to low-frequency excitations while (b) and (d) correspond to high-frequency 
excitation. 
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FIG. 6. Nonlinear Kerr rotation angles for p polarized incident light (p^Kp (full and short-dashed 
curves) and for s polarized incident light (long-dashed and dotted curves) for Fe at 770 nm 

as a function of the angle of incidence 6 in the longitudinal Kerr configuration. The relative phase 

(2) (2) (2) 

between Xzxx = Xzyy and Xzzz is </5i=0.5057r in the full and long-dashed curves and <^2=l-5057r 
in the short-dashed and dotted curves. 
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